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Abstract. We describe an operational scheme for determining both the position and momen- 
tum distributions in a large class of quantum states, together with an experimental implemen- 
tation. 
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! The basic understanding of quantum mechanics asserts that due to their noncommutativity, 

$_i | position and momentum of a quantum object cannot be measured together in any state. It is 
an equally well known fact that these observables can be approximately measured together if 
the measurement accuracies fulfill an appropriate trade-off relation. Much effort has gone into 
obtaining rigorous formulations of these facts; for a recent overview see, for instance, pp. 

In this paper we wish to point out that there are measurement schemes which allow one to 
determine both the position and the momentum distributions from the statistics of a single 
measurement: in this sense the position and momentum of a quantum object can, after all, be 
measured together. 

Technically, this observation builds on the well-known properties of some covariant phase 
space observables, such as the Husimi Q-function. The measurement in question can then be 
implemented via eight-port homodyne detection. Therefore, the idea is completely realistic. 

The structure of covariant phase space observables is well-known. If S is a positive operator 
of trace one (acting on the Hilbert space TC of the quantum object), and W qp , (q,p) G M 2 , are 
the Weyl operators associated with the position and momentum operators Q and P, then the 
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(with the integral being defined in the sense of the weak operator topology) constitute a covari- 
ant phase space observable G s (defined on the Borel subsets of the phase space IR 2 ). Conversely, 

£a! \ any covariant phase space observable is of that form for some S [2]. If the generating opera- 
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tor 5* is the vacuum state |0), then the probability distribution Z i— > G -AZ) = Tr[pG (Z)] 



Tr[pG'°^(Z)] is the Husimi distribution of the state (density operator) p. 



We are going to show that for a certain family of operators S the phase space observable 
G s entails complete information about the position and momentum observables: that is, the 
spectral measures of Q and P can be determined from G s . This is done by using the moment 
operators L(x k , G s ) and L(y k , G s ) of G s , which are defined as follows (using the general theory 
of operator integrals) : 

L(x\G s ) = -L / q k W qp SW*dqdp, 

L{y\G s ) = f p k W qp SW* qp dqdp. 

In general, it is very difficult to determine the domains of these (unbounded) operators. How- 
ever, the following is known [?l Theorem 4] (for a similar result, see |8j): for any k G N, if (and 
only if) the operator Q k \fS is a Hilbert-Schmidt operator, that is, Tr[5 , Q 2fc ] < oo, then 

k 

(1) L(x k ,G s ) = J24Q\ with si = (*)(-l)*-^[g*-'5], 

1=0 



and, similarly, if (and only if) the operator P k \fS is a Hilbert-Schmidt operator, then 

k 

(2) L(y k ,G s ) = J2 s M pl > with s ki = 0(-l) k - l Tr[P k - l S]. 

1=0 

For any vector state <p, the numbers (ip\L(x k , G s )<p) are the k th moments of the marginal 
probability measure X \— > (ip\G s (X x R)</?), and similarly for the numbers (<p\L(y k ,G s )(p). 
From Eq (PQ) one can easily determine any Q k in terms of the (operationally well defined) 
moment operators L(x l , G s ), I — 0, . . . k, and, similarly, each P k in terms of the (operationally 
well defined) moment operators L(y l , G s ), I = 0, . . . k, from ([2j). More explicitly, consider the 
simple recursive transformation 

fc-i 

(3) Pk = a k -^s kl (3 h ken. 

1=0 

Putting a k = (f\L(x k , G s )<p) and s k i = yields [3 k = ((p\Q k <p). (Note that s® k = 1.) Similarly, 
one gets the moments (ip\P k <p) from the actual marginal moments ((p\L(y k ,G s )(p). Therefore, 
the numbers (ip\Q k ip) and (ip\P k <p) are just as operational as the actually measured moments 
(ip\L(x k , G s )ip) and {ip\L(y k , G s )(p), from which they can be reconstructed via ([31). But (ip\Q k ip) 
and (ip\P k ip) are the moments of the probability distributions p® and of the position and 
momentum observables Q and P in the initial state <p. If (p is chosen to be, for example, a 
linear combination of Hermite functions, the measures p® and p^ are exponentially bounded, 
and, as such, are uniquely determined by their respective moment sequences {{<p\Q <p))k an d 
(((p\P k ip)) k p. 406, Theorem 30.1]. In this sense one is able to measure the noncommuting 
observables Q and P simultaneously in such a vector state (p. Furthermore, since the linear 
combinations of Hermite functions are dense in £ 2 (R), their associated distributions p® and 
p^ suffice to determine the whole position and momentum observables Q and P (as spectral 
measures). 

The phase space observables G s , with the generating operator S being such that Tr[W gp S] ^ 
for almost all (q,p) G M 2 , are informationally complete: that is, the statistics G s p determine 
the state p of the quantum object. Given the state p, the distributions of all observables are 
then determined, including, in particular, p® and p p p . In this sense one might argue that an 
informationally complete measurement serves as a measurement of any other observable of the 
system. However, to our knowledge, there is no simple method to determine from the statistics 
G p the statistics of another observable. An explicit reconstruction of the state typically requires 
several different measurements, as is the case with the usual tomographic methods. 

In the present case we are able to construct directly and in simple operational terms (without 
any use of the informational completeness) the moments of the distributions^ and p p from the 
statistics G^ for any state ip, and these moment sequences identify, in turn, the distributions 
themselves in the exponentially bounded case. 

It is useful to note that the moment sequence of, say, the momentum distribution does not, in 
general, identify that distribution, even in cases where all the moments are finite. For example, if 
ips = (ipi + e lS ip2) / V2 is a double-slit state, that is, ^1,^2 are (smooth) wave functions supported 
on disjoint intervals, then the moments ((ps\P k <ps) of the momentum distribution p^ s are all 
finite and independent of 5 although the distribution does depend on 5. This highlights the 
fact that the assumption of exponential boundedness is crucial in the above consideration. 

We also want to point out that the exponential boundedness condition itself is operational, 
since it can, in principle, be verified by considering the very moment sequence that is obtained 
in the measurement. In fact, a probability measure with the moment sequence m 1 ,m 2 , ... is 
exponentially bounded if and only if 

(4) \m k \ < CR k k\, k = 1,2,..., 

for some C, R > (see e.g. the proof of Proposition 2 of |10j). 
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Experimental implementations of observables G can be obtained in quantum optics. For 
example, the Husimi Q-function of a single-mode signal field can actually be measured by 
the eight-port homodyne detection technique, provided that the reference beam (used for the 
homodyne detection) is a very strong coherent field so that it can be treated classically; see 
e.g. [H]. In a recent paper it has been shown, without any classicality assumptions, that any 
covariant phase space observable G s can be realized as a high-amplitude limit of an eight-port 
homodyne detection observable choosing the auxiliary parameter field to be in the conjugated 
state C^SC 02]. 

In this measurement, the input signal field, in state p, is mixed by a beam splitter with another 
single-mode field, the latter being in the reference state C~ 1 SC. Each of the two output beams 
of the beam splitter are led to a homodyne detector which measures a field quadrature. One 
measures -^(a* + a) and the other one ^i(b* — b), where a and b are the annihilation operators 
for the signal and the reference field (in the Schrodinger picture). When this measurement is 
interpreted with respect to the signal field in the initial state p, the combined statistics of the 
homodyne detectors are given by the distribution (X,Y) t— > Tr[pG s (X x Y)], where the Weyl 
operators W qp are now defined in terms of the signal field quadratures Q := ^(a* + a) and 

P:=^i(o*-o). 

By performing the measurement, one thus obtains two outcomes, which then are to be 
interpreted with respect to the signal beam. Notice that the outcomes for the two homodyne 
detectors are obtained at the same time. Consequently, upon repeating the preparation and 
the measurement, the statistics for the detectors are likewise created simultaneously. After a 
suitable statistical accuracy has been obtained, and the measurement is thus complete, one 
takes the statistics given by, say, detector 1, computes its moments a^, and calculates the 
numbers /3& using ([3]). The f3 k are now the moments of the distribution of Q, in the input 
state of the signal field. A similar procedure applied to detector 2 yields the moments j3' k of 
the distribution of P. One can then check whether (pj holds for both the sequences (3k and 
(3' k . If that is the case, then it is justified to say that (the statistics of) the position- and 
momentum-like observables Q and P have been measured together. 
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